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^ ■ Abstract 

I We present an analysis of QCD sum rules for pion form factor in next-to-leading order of 

perturbation theory for the case of axial- vector pion currents. The theoretical predictions for 
C^c Q 2 -dependence of pion form factor are in good agreement with available experimental data. 

It is shown, that NLO corrections are large in this case and should be taken into account in 
£NJ ' rigorous theoretical analysis. 

(N ' 

O ! 1 Introduction 

■ The study of electromagnetic form factors of hadrons already has a long history. The first 

applications appeared right after the realization, that perturbative QCD may be also applied 
in studies of exclusive processes with high momentum transfer P3E3IH1- However, later compar- 
ison of pQCD predictions with data lead to conclusion, that at moderate momentum transfers 



(typically of the order of few GeV) the power corrections, otherwise known as soft or end-point 
contributions should come into play 1 . 

In this paper we are using the framework of QCD sum rules jl] to consistently account 
for both hard scattering and soft wave function overlap contributions to pion electromagnetic 
form factor. Within this approach the soft contribution is dual to the lowest-order triangle 
diagram, while the hard contribution is given by diagrams having higher order in coupling 
constant a s and as a consequence suppressed relative to soft contribution with additional factor 
a s /ir ~ 0.1. This extra suppression is in a complete agreement with asymptotic behavior of the 
pion electromagnetic form factor, calculated in pQCD ^|21|3]: 
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where the last equality holds for the asymptotic pion distribution amplitude. At asymptotically 
high Q 2 the O(o s /vr) suppression of hard contribution is more than compensated by its slower 
decrease with Q 2 . However, such compensation does not occur in the region of moderate 
momentum transfer, where the soft contribution, scaling as l/Q 4 becomes large and can compete 
in strength with hard contribution. 

The pion electromagnetic form factor was studied using many different frameworks, like QCD 
sum rules |BJ El El E] , light-cone sum rules ^3 ^JJ , sum rules with nonlocal condensates [T5] 
and NLO pQCD approach [141 1151 [TBI 1171 118j . improved by inclusion of finite size corrections 



1 See, for example, discussion in 
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to pion distribution function |19l 1201 I21| . There are also estimates of pion electromagnetic form 
factor based on use of pseudoscalar pion interpolating currents [2H HH1 HI] • 

In what follows we will consider NLO three-point QCD sum rules for pion form factor, where 
axial currents are used as pion interpolating currents. The main result of this paper is the explicit 
analytical expression for QCD radiative corrections to double spectral density, entering QCD 
sum rule predictions for pion electromagnetic form factor. It should be noted, that up to the 
moment there are only radiative corrections computed for reduced spectral density |25[ I26j . 

The obtained results for pion electromagnetic form factor calculation are in good agreement 
with existing experimental data. Here, we would like to stress that an inclusion of radiative 
corrections is very important, as only in such a way we can simultaneously account for both 
hard and soft contributions. Moreover, these corrections are large numerically and thus should 
be accounted for in theoretical predictions. 

The paper is organized as follows. In section 2 we describe our framework and give explicit 
expressions for next-to-leading order corrections to double spectral density. Section 3 contains 
our numerical analysis. Finally, in section 4 we draw our conclusions. 

2 Derivation of QCD sum rules 

To determine pion electromagnetic form factor we will use the method of three-point QCD sum 
rules. Here, to describe charged pion state we choose an action of axial interpolating current on 
a vacuum state. The vacuum to pion transition matrix element of axial current is defined by 

(0|n757^d|vr _ (p)) = if^, (2) 

where f n = 131 MeV. Next, the pion electromagnetic form factor, we are planning to determine, 
is given by hadronic matrix element of electromagnetic current: 

WM 1 Wp))=fAQ 2 )(p»+p'^ (3) 

where j^ 1 = e u u^^u + edd^^d, the momenta of initial and final state pions were denoted by p, 
p' and Q 2 = —q 2 (q = p — p') is square of momentum transfer. 

Within the framework of QCD sum rules an expression for pion electromagnetic form factor 
follows from an analysis of corresponding three-point correlation function: 

U fmP {p,p' : q)=i 2 J ^ e i ^--P^(0|r{n(x)757 Q d(^),^ I (0),(n(y) 7 57^(y)) + |0) (4) 




Figure 1: LO diagram 



This correlation function contains a lot of different tensor structures. The scalar amplitudes 
Hj in front of different Lorentz structures are the functions of kinematical invariants, i.e. = 



2 



Ui(p 2 ,p' 2 , q 2 ). The calculation of QCD expression for three-point correlator is done through the 
use of operator product expansion (OPE) for the T-ordered product of currents. As a result 
of OPE one obtains besides leading perturbative contribution also power corrections, given by 
vacuum QCD condensates. We will return to the discussion of QCD expression for three-point 
correlation function in a moment. Now let us discuss the physical part of our sum rules. The 
connection to hadrons in the framework of QCD sum rules is obtained by matching the resulting 
QCD expressions for current correlators with spectral representation, following the structure of 
double dispersion relation at q 2 < 0: 



rr I 2 2 2\ 
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(2ir) 2 J ( Sl -pl)(s 2 -pl) 



ds\ds2 + subtractions. 



(5) 



Assuming that the dispersion relation © is well convergent, the physical spectral functions 
are generally saturated by the lowest lying hadronic states plus a continuum starting at some 
thresholds sf 1 and s 2 h : 



9{si _ s f ) . 9{S2 _ s f) . p conJ (si)S2)Q2)) 



(6) 



where 



(0\u l5la d\ IT' (p')) <7T- {p') | jl V (P)> <T~ (P) I (n757/3^) + 1 0) 

(2tt) 2 S(si)5(s2) + higher state contributions 



(7) 



In our approximation of massless quarks we put m 2 = 0. So we see, that pion contribution 
into spectral density is given by p^°p ~ f 2 F 7T (Q 2 )p a p' l3 (p fl + p ,fl ) and as was already noted 
in [HJ |S1 EI the most convenient way to extract pion form factor from QCD sum rules is to 
consider scalar amplitude in front of most symmetric Lorentz structure P^P a Pp {P =p + p'). 






Figure 2: NLO diagrams 



Now it is time to return to the problem of evaluation of QCD expression for three point 
correlation function, we are interested in. The condensate contribution is known already for a 
long time jSJ El E] and its analytical expression could be found in the section with numerical 
results. The calculation of perturbative contribution could be conveniently performed with the 
use of double dispersion representation in variables s\ = p 2 and s 2 = p' 2 at q 2 < 0: 

npert, 2/2 2^ 1 f P^( S 1^2,Q 2 ) 

n£ a/3 (p ,p ,q) = J ( Sl _ p 2 )(52 _ p /2) & i & 2 + subtractions (8) 
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The integration region in (jSJ) is determined by condition 2 

- 1 < ^~ Sl ~ g2 < 1 (9) 

and 

\(XI, X 2 , X 3 ) = (Xl +X 2 - X 3 ) 2 - 4X1X2. (10) 

The double spectral density P^p(s\, S2, Q 2 ) is searched in the form of expansion in strong 
coupling constant: 

*2, Q 2 ) = P%(si,s 2 , Q 2 ) + (g) , 2 , Q 2 ) + . . . (11) 

At leading order in coupling constant we have only one diagram depicted in Fig. 1, con- 
tributing to three-point correlation function. At next to leading order we have 6 diagrams shown 
in Fig. 2. The calculation of corresponding double spectral density was performed with the 
standard use of Cutkosky rules. In the kinematic region q 2 < 0, we are interested in, there is no 
problem in applying Cutkosky rules for determination of p^ap(si, S2, Q 2 ) and integration limits 
in s\ and S2- The non-Landau type singularities, not accounted for by Cutkosky prescription, 
do not show up here. 

It is easy to find, that at Born level the scalar spectral density in front of most symmetric 
Lorentz structure P^P a Pp is given by: 

A(S1,«2,Q 2 ) = ^^(mSl + S2 + Q 2 )(si+S 2 + 2Q 2 )-k 2 



5Q 2 (s 1 + s 2 + Q 2 ) 3 )p ix P a P l3 + (12) 



where k = A(si, S2, —Q 2 )- The full analytical expression for p^p could be found in 0. The cal- 
culation of NLO radiative corrections to double spectral density is in principle straightforward. 
One just needs to consider all possible double cuts of diagrams, shown in Fig. 2. However, the 
presence of collinear and soft infrared divergences calls for appropriate regularization of aris- 
ing divergences at intermediate steps of calculation and makes the whole analytical calculation 
quite involved. We will present the details of NLO calculation in one of our future publications. 
Here we give only final results. The calculation could be considerably simplified with the help 
of Lorentz decomposition of double spectral density based on a fact, that our spectral density 
is subject to three transversality conditions: p^a^Qfi = P^apPa = P^app'p = ; 

pH°P = A, [(Q 2 + xK + y)j%\ \{y - x)j4 + (Q 2 + x)p$] [{Q 2 + y)ft + (Q 2 - y)p 2 l ] 
-\a 2 [(Q 2 + yK + (Q 2 - yK][(Q 2 + x)g<# - 2p^} 

-\a 3 [(q 2 + xK - (x + vMMpI - pfK + (Q 2 + V)9^] 

-\m[{x - y)p{ - (Q 2 + x)p%\[2(j% - + (y - QV«], (13) 

where x = s\ + S2 and y = s\ — s 2 - The four independent structures A4 (we suppressed the 
dependence on kinematical invariants) are given by a solution of system of linear equations: 

h = p m pp1p^P 1 =^(kA 1 -A 2 -A 3 -AA (14) 



2 In our case this inequality is satisfied identically. 
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h 



\{x + Q 2 ) 



kAi - ?>A 2 - A 3 



P^a/3Pl9^ 



k 



(y + Q 2 )[kA 1 -A 2 -3A 3 -A 4 



= --(y- Q 2 ) [ kA l -A 2 -A 3 -3A 4 ), 



(15) 
(16) 
(17) 



The analytical expressions for Ii (functional dependence on kinematical invariants is assumed) 
were found to be (53 = Q 2 ): 



k l l 2 h 



-s\ + S 2 s\ + slsi - S% + (Si + S 2 )sl - s\s 3 - s\s 3 + SlS 2 S 3 

-161og 2 («i) - 161og(w 3 )log(t>i) - 161og(u 4 )log(vi) 

+21og(vi) - 41og 2 (v 3 ) - 41og 2 (w 4 ) - 21og(w2) - 21og(v 3 ) - 81og(v 3 ) log(« 4 ) 



-8Li 2 



8Li 2 



8Li 2 



8Li 2 



+ 8Li 2 



k l ' 2 h 



-2s\ - 2s\ + 2s\ - 8sis 2 + s\s 2 



-32 log 2 M - 321og(v 3 )log(vi) - 321og(v 4 )log(«i) + 41og(m) 
-81og 2 (v 3 ) -81og 2 (t> 4 ) -41og(u 2 ) -41og(u 3 ) - 161og(v 3 )(v 4 ) 



- lliLi ? ( — 



i.()Li, ( y± 

V2 



16Li 2 I ^ 



si 



16Li 2 ( — ) + 16Li 2 ( — 

Si 



k l ' 2 h 



-2s\ + 2s 2 2 + 2s 2 - 8s 2 s 3 + s 2 s 3 



-321og 2 (t;i) -321og(w 3 )log(ui) - 321og(« 4 )log(vi) +41og(«i) 
-81og 2 (v 3 ) - 81og 2 (u 4 ) - 41og(v 2 ) - 41og(u 3 ) - 161og(t> 3 ) log(v 4 ) 



- KiLi , ( — 

Xl 



KiLi, [ y± 

V2 



KiLi , I — 

si 



IGLi, ( £1) +l6Li 2 ( - 

si) \z 2 



k l l 2 I A 



2s\ — 2s\ + 2s 3 — 8S1S2 + s\s 3 



-32 log 2 M - 321og(?; 3 ) logfa) - 321og(« 4 ) logfa) + 41og(m) 
-81og 2 (w 3 ) - 81og 2 (t> 4 ) - 41og(t> 2 ) - 41og(u 3 ) - 161og(t> 3 ) log(u 4 ) 



-16Li 2 



16Li 2 



16Li 2 



- 16Li 2 



16Li 2 



where the following notation was introduced: 



(18) 



(19) 



(20) 



(21) 



Xl = 


2^1 - si- Q 2 ) - ^Vk, 


(22) 


x 2 = 


\{si-s 2 -Q 2 ) + -Vk, 


(23) 


yi = 


^{si + Q 2 - si) - -Vk, 


(24) 


V2 = 


^( Sl + Q 2 -s 2 ) + ^Vk, 


(25) 


Z\ = 


^(si + s 2 + Q 2 )-^Vk, 


(26) 


z 2 = 


^(si + s 2 + Q 2 ) + -Vk, 


(27) 


Vl = 




(28) 



t>2 = 


1 

2s 2 


(si -s 2 + Q 2 ) + 


2s 2 


(29) 


7 Jo 


1 




— 


f30) 

V ou 7 




s\ 
Q 2 ' 






(311 


^5 = 


Q 2 ' 






(32) 


V 6 = 


1 - 


21 




(33) 
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We checked, that all infrared and ultraviolet divergences cancel as should be for axial interpolat- 
ing currents. Finally, NLO scalar spectral density in front of most symmetric Lorentz structure 

P^ip = ^l^( x i + x 2)kl3 + (k-5(x 1 +x 2 )(yi+y 2 ))Ii 

Ilk 1 

+^(yi + V2)kh + - — ■ ((xi + x 2 )(yi + y 2 ) - k)I 2 > P^P a Pp + ... (34) 

2 2 zi + z 2 J 

In the limit Q 2 — ► 00 our NLO double spectral density takes the following form: 

Pile = {q2- w ^t~ ~ i^r- log J j W„ + • ■ ■ (35) 

Here we would like to make several comments. First, we see that double logarithms are absent in 
our answer. Typically we would expect them to be nonzero as diagrams in FigEJcontain Sudakov 
vertex - corrections to g-vertex. And indeed our result for gluon correction to electromagnetic 
vertex (accompanied by the appropriate 1/2 self-energy insertions to ensure ultraviolet-finite 
result) contains double logs and agrees with large Q 2 limit of |25| I26j . However, our results for 
corrections to p\ and p 2 vertexes also contain double logs (an analogue of double logs arising 
in pQCD description of pion electromagnetic form factor |19l !2L)j ) . In the sum of all diagrams 
double logs cancel. Second, using large Q 2 expression for our spectral density and continuum 
subtraction thresholds for s\ and s 2 equal to 47r 2 / 2 it is easy to see that leading Q 2 — » 00 QCD 
sum rule expression for pion electromagnetic form factor is given by leading pQCD prediction 
(^Q) with asymptotic pion distribution function (c/>" s (x) = 6/ vr x(l — x)) . Third, we checked 
numerically, that a Q 2 — > limit of our spectral density do satisfy Ward identity constrains 3 
|27| . Now let us proceed with numerical analysis. 



3 Numerical analysis 



In numerical analysis we used Borel scheme of QCD sum rules. That is, to get rid of unknown 
subtraction terms in (JSJ) we perform Borel transformation procedure in two variables s% and s 2 . 
The Borel transform of three-point function Hi(s\, s 2 ,q 2 ) is defined as 



<S>(Mt,M 2 ,q 2 ) = B 12 U i (s 1 ,s 2 ,q 2 " 



sl +1 ( d\ n s™ +1 f d 



n,m— »oo ) //! \ ds 2 J ml V dsi 



si=mMj ! ,s 2 =nM% 



Tl(s 1 ,s 2 ,q 



2\ 



(36) 



3 We are grateful to A.P.Bakulev for explaining to us the details of both this limit and Ward identity constrains for 
double spectral density 
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Then Borel transformation (|36|) of (jHJ) and © gives 



$ (pert|phys)^ M 2 )M 2 )g 2 



(2tt)s 



ds 2 exp 



si 

M 2 



■S2 



/9 (pGrt|phyS) (si,«2,9- 



2 \(37) 



where p( pert \ phys \si, S2,q 2 ) stands for the expression of scalar spectral density in front of 
P^PaPp Lorentz structure. In what follows we put M 2 = Mf = M 2 . If M 2 is chosen to 
be of order 1 GeV 2 , then the right hand side of (|37|) in the case of physical spectral density will 
be dominated by the lowest hadronic state contribution, while the higher state contribution will 
be suppressed. 



FJ1 GeV 2 
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Figure 3: Borel mass M 2 dependence of pion electromagnetic form factor at Q 2 = 1 GeV 



Equating Borel transformed theoretical and physical parts of QCD sum rules we get 
FAQ 2 ) = j, (*(MV) + ^(OI^G^IO) + ^a s (0|^|0) 2 (l + ^ )) , (38) 
where 

f X dypV ert ( Sl , S2 ,q 2 ). (39) 


Here, for continuum subtraction we used so called "triangle" model. To verify the stability of 
our results with respect to the choice of continuum model we checked, that the usual "square" 
model gives similar predictions for pion electromagnetic form factor provided so ~ 1.5si is 
chosen 4 . In what follows we use sq = 0.96 GeV 2 for continuum threshold 5 . This value is in 
agreement with the continuum threshold ~ 0.7 GeV 2 for axial polarization operator used in 
two-point sum rules. Next, we use two-loop renormalization group running of strong coupling 
constant with Aqcd = 325 MeV and fix the scale \x of strong coupling constant at 2 GeV. This 
choice is in agreement with the discussion presented in |2H], where it was argued that in the 
region of momentum transfers Q 2 < 10 GeV 2 the strong coupling constant should be taken 

4 For more information about different continuum subtraction models see [2] 

5 In general, the value of continuum threshold is determined from the ratio of nonperturbative contribution to 
leading perturbative term in OPE for our correlation function. 



$(M 2 ,g 2 ) 



dx exp 



x 
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at frozen value a s ~ 0.3. In Fig. |3] we plotted the dependence of the pion electromagnetic 
form factor from the value of Borel parameter M 2 at Q 2 = 1 GeV 2 . It is seen that "stability 
plateau" starts to develop for M 2 > 2 GeV 2 . To proceed further we could fix the value of Borel 
parameter at M 2 = 2 GeV 2 and get results for pion form factor as a function of momentum 
transfer. However, doing so we limit the range of Q 2 where our results could be considered as 
reliable. To see it, it is instructive to investigate large Q 2 limit of pion form factor. In the 
limit Q 2 — ► oo perturbative contribution to pion form factor decreases as Q~ 2 , while power 
corrections grow with Q 2 . It turns out, that our sum rules become inapplicable at relatively 
large momentum transfers Q 2 > 4. Therefore, to cover all experimentally accessible range of 
Q 2 for pion electromagnetic form factor we will explore our sum rules in the limit of infinite 
Borel parameter M 2 . So, basically, here we are employing local duality approach. The local 
duality approach implies the following relation for continuum threshold |2fi[ 127): s\ = S2 = 
4tt 2 / 2 /(1 + St) = 0.62 GeV 2 , which ensures the Ward identity for pion electromagnetic form 
factor up to NLO (i^(0) = 1). The results for pion electromagnetic form factor are shown in 
Fig. 21 (solid line is the sum of LO and NLO contributions, curve with long dashes denotes LO 
contribution and curve with short dashes stands for NLO contributions.). 




Now, let us recall that in the limit Q 2 —* oo our QCD sum rule prediction coincides with 
LO pQCD results with asymptotic expression for pion distribution function. We could improve 
our prediction further with inclusion of NLO pQCD corrections. 

Within perturbative QCD pion electromagnetic form factor is given by the following factor- 
ized expression 



ppQCG 



dx I dyfaix, n)T H (x,y,Q 2 , fj?)4>n(y, fi) 



(40) 



with 



TH{x,y,Q 2 ,n 2 ^ 



27rC F a s ^)f 2 
N C Q 2 (1 - x)(l - y) 



1 



a s (n) 



T 1 (x,y,Q 2 / f i 2 ) + 



IT 



(41) 



It should be noted that this contribution contains explicit dependence on pion wave function and 
gives us a possibility for its determination through the comparison between combined predictions 
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Figure 5: Q 2 dependence of pion electromagnetic form factor including NLO pQCD corrections 



of QCD sum rules and pQCD frameworks with experimental data. We suppose to perform this 
analysis in future, while in present work we limit ourselves only to the inclusion of NLO pQCD 
prediction with the use of asymptotic pion distribution functions. Within this approximation 
pQCD predictions for pion electromagnetic form factor are given by |14[ I15| I16| I17[ I18j : 

F pQCD (Q2) = F LO (g2) +i? NLO (Q 2 ); 



where 



F LO (Q 2 ) 
F NLO (Q 2 ) 



f 2 



2 

Q 2 ' 



6.58 + - log 



Q 2 



(42) 

(43) 
(44) 



The results for pion electromagnetic form factor including NLO pQCD corrections are shown in 
Fig. 0] (solid line is the sum of QCD sum rule and NLO pQCD contributions, curve with long 
dashes denotes LO QCD sum rule contribution and curve with short dashes stands for the sum 
of NLO QCD sum rule and NLO pQCD contributions.). We see, that obtained results for pion 
electromagnetic form factor are in good agreement with available experimental data as well as 
already available theoretical estimates performed with the same level of precision. 



4 Conclusion 

We presented the results for pion electromagnetic form factor in the framework of three-point 
NLO QCD sum rules with axial interpolating currents for pions. To improve our estimates we 
also took into account NLO pQCD contributions with asymptotic pion distribution functions. 
The theoretical curve obtained for Q 2 dependence of pion form factor is in a good agreement 
with existing experimental data. Here, we for the first time computed radiative corrections to 
double spectral density entering three-point sum rules with axial currents. These corrections 
turned out to be large and should be taken into account in rigorous analysis. 
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